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Abstract. Planar solidification from an undercooled melt has been considered 
using the phase-field model. The solute and the phase fields have been found in 
the limit of small impurity concentration. These solutions in the limit of vanishing 
velocity of the interface motion give the equilibrium partition coefficient and the 
liquidus slope. Asymptotic expansions for the solute and for the phase fields, and 
the relation between the diffusive speed and the parameters of the phase field model 
have been found at high growth velocity. A comparison with numerical calculations 
is presented. 

1 Introduction 

Classical macroscopic models of solidification are based on the assumption 
that the interface between regions of different phases have zero thickness pQ . 
Sharp-interface descriptions require the introduction of separately derived 
nonequilibrium models for the behavior of the interfacial temperature and 
of the solute concentrations. In contrast to the sharp-interface models, the 
phase-field models describe the bulk phases as well as the interface, i.e. 
these models treat the system as a whole and eliminate the need to specify the 
interfacial conditions separately. In phase-fields models, equilibrium behav- 
ior is recovered at low growth velocities and nonequilibrium effects naturally 
emerge at high growth velocities. From a computational point of view, the 
advantage of the phase-field formulation is that the interface is not tracked 
but is given implicitly by the value of the phase-field variable. However, the- 
oretical analysis solutions were obtained by matching of separately derived 
"inner" (near the interface) and "outer" (far from the interface) solutions 
in an intermediate region 1 1I5| . In the work phase-fields and solute pro- 
files have been obtained without using the multiply-variable expansion, which 
leads to "inner" and "outer" solutions, but the profiles were obtained from a 
simplified set of equations in which the interface kinetics was eliminated. 

In this work an analytical solution of the phase-field model for planar 
solidification from an undercooled melt as an expansion in terms of the solute 
concentration far from the interface is presented. This solution is valid for the 
bulk phases as well as for the interface and includes the interface kinetics. To 
verify the expansion, numerical simulations of the problem were carried out. 
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2 The Model 

Let us consider the motion of a planar solid-liquid interface during the solid- 
ification of a binary alloy. The governing equations for the phase-field <fi(x, t) 
and for the solute concentration c[x, t) are given by the following set of equa- 
tions [2] 

d<p AT ( 2 d 2 <t> df\ dc a d df 

which satisfy the conservation of solute and which have been derived under 
the assumption that the total free-energy decreases monotonically in time. 
The parameters Mi, M2, and e are related to the growth kinetics, to the 
diffusion coefficient, and to the surface energy respectively. 

In the approach of ideal solution, the free-energy density /(</>, c, T) at a 
temperature T is given by [2] 

f{d>, c, T) = cf B (4>, T) + (1 - c)f A (4>, T) + — [clnc+ (1 -c) ln(l - c)], (2) 

where R is the universal gas constant, v m is the molar volume, which is 
assumed to be constant. Free-energy densities of the pure materials A and B 
are assumed to be of the form |2] 

fifaT) = Wi /"%(?-!) (p-l- Pi(T)) dp = ^g^+^i^lp^) (3) 



where 



5 (0)=0 4 -20 3 + 2 , p(0) = ^(3-20), ^fl = Li ^- *L. (4) 

6 Li 

Here Wi is a constant, Li is the latent heat per unit volume, Tj is the melting 
point of material i = A, B. 

To investigate steady-state interface motion, the boundary conditions will 
be taken in the form 

dc 

C \x — >oo ^oc, $\x — >oo 0, ~7\ \x — > — oo 0, 4*\x — > — 00 1, (5) 

OX 

and the moving frame z will be adopted according to the equation 

z = x - Vt, (6) 



so that the interface (given by = 1/2) corresponds to z = at this frame. 
In the moving frame, equations Q may be written by using as 
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M 2 RT d 2 c 

v m dz 2 



dz 



M 2 — 

dz 



C{ C> 1 ^ dcj) J dz 



= 0. 



(8) 



For further analysis we introduce the notation for the diffusion coefficient 
D = MiRT jv m and integrate equation JHJl with respect to z over an interval 
[z, oo). Employing the boundary conditions |pj gives 



Dp + V(c - Cco ) + M 2 c(l 



dj)_ 
~d~z 



0. 



(9) 



3 An Approximation of a Small Impurity 
Concentration 

Let us examine a dilute binary alloy with initial impurity concentration <C 
1. For this purpose we shall consider the expansion in terms for the phase- 
field, the solute concentration, and the temperature 

<j>{z) = <t> {z) + c oo 4> 1 {z) + O{(? oo ), (10) 
c(z)= CooCl (z) + 0(cL), (11) 
T = T + CocTx + 0(4), (12) 

As a result, in the zeroth-order approximation, from by using J3J we 
obtain the equation for (j>o(z) 

V d<t> {Wa£>i 
6 dz 2 Mi dz I 4 <90 

The solution of (|13fl . which satisfies the boundary conditions (0, describes 
the phase field during the solidification of a pure material A and has the form 



(b<\ 



W a (3a dp 
6 dd> 



= 0, 



(13) 



M z ) = 2 



1 — tanh 



-) 

21 A J 



(14) 



where I a is the interface thickness and the interface velocity V is related to 
the temperature T by the relationship 



V 



6MiIaLj 
Ta~ 



(Ta-To), 



(15) 



which determines the interface kinetics for a pure material. The relationship 
1151) corresponds to the model of collision limited growth [JJ with a kinetic 
coefficient 



/' 



QM^aLa 



T = T a 



V 



(16) 
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In the first-order approximation, the equations (Q) and © give 



v dfa 

dz 2 Mi dz 



d 2 f A 



4>o,Ta 



Cl 



rp d 2 f A 






1 dTd<p 


4>o,T 




(df B df A \ 




V d<t> 


dcj> ) 


00, T 



0, (17) 



dci V, OAF n 

where function Zii^z) is defined by expression 



(18) 



AF{z) 



RT 



(19) 



Taking into account the boundary conditions JHJ), the solution of ljl8|l is 

Cl ( z ) = ^ T e ^(«'-»)/o+^'(*')-^'W dz ' 



(20) 



The point z = corresponds to the solid-liquid interface, i. e. 0(0) = 
^o(0) = 1/2, and to the inflection of the phase-field, therefore 



n\z=o 



0, 



dz 2 



= 0. 



(21) 



As for metallic systems Mi ~ 10 8 cm 3 /J-s [2], then the relation V/M\ <C 1 
will be satisfied. Taking into account the conditions l|21(l at z = 0, one can 
find from H17fl the relationship between T\ and Ci 



Tx = - 



(df B df A 

\ dcj> d(j> 



d 2 f A 




dTd4> 


<t>o,Ta 



ciU=o- 



(22) 



Equations fTTft . (ffT^f) . and (|2^|l give the expression for the temperature 



T = T A m(V)c a=0 , 

A* 

where the function m(V) is defined by 



(23) 



m(V) 



W B /3 B (T)-W Ay 9 A (T) 







dTd(/> 


4>o, T 



w A 



d(3 A 
dT 



(24) 



T=T A -V/n 
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4 Stationary Interface 



Now we consider the solute field due to the stationary solid-liquid interface, 
at V = 0, which corresponds to the equilibrium state. Let us introduce the 
definitions for the bulk concentrations in the liquid cl = c(z — > oo) and in 
the solid c$ = c(z — > — oo). For a stationary interface, first it follows from 
l(TB|) that T =Ta, and second the equation gives the solute field 



£(£) 

CL 



which leads to the equilibrium partition coefficient 



(-L 



exp[— AF(z — > — oo)] = exp I — 



v m L B {TA - T B ) 



RTaTb 



(25) 



(26) 



This expression has the same form as in [H], where the equilibrium partition 
coefficient has been obtained from a simplified set of equations in which the 
effects of interface attachment kinetics are eliminated. 
Using Q), from equation l|24|l at To = Ta we obtain 



too = m(0) = j^t^( t a - T B ). 



(27) 



The solute concentration at the interface z = 0, in accordance with (|25|l . is 

' W b -Wa , 1L B (T A -T B ) 



c\z=o = c L exp 



RT A 



64 



(28) 



Thus equations l)23[l. I|27|) . and l|28|) give the temperature at the equilibrium 
solid-liquid interface 



T = T A - m c| z=0 = Ta - m exp 



+ 



IL B (T A -T B ) 



Wb-Wa 
64 



c L =T A - m e c Ll (29) 



2 T B 

which yields the following equation for equilibrium liquidus slope 



Lb Ta , , 
— — {T A -T B )exp 



RT A 



Wb-Wa , L B (T A -T B ) 



64 



2T B 



(30) 



5 Large- Velocity Asymptotics 



To investigate the behavior of the solute concentration under rapid solidifica- 
tion conditions, we consider an asymptotic analysis in the limit F>1. The 
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solution 120(1 can be written in the form 

Cl ( z ) = L [°° e -Vz'/D+AF {z - z ')-AF( z ) dz ^ 

D Jo 



(31) 



Applying the Laplace's theorem about asymptotic expansions [8] to this in- 
tegral, we obtain the solute field with an accuracy to second order 



d(z) 



c(z) 



1 



DdAF 
V dz ' 



(32) 



According to 1)32(1 . a maximum value of solute is reached at z = 0. Therefore 
the segregation coefficient at large velocities is given by 



k = 



1 



C\z=0 



D dAF 



V dz 



D dAF 

= 1 + — 



z=0 



z=0 
= 1 - 



V dz 
D v m ( df B df A 



Al A V RT 



(33) 



i=l/2 



For the sharp-interface model the nonequilibrium segregation coefficient 
has been derived by Aziz [5] and in the limit of dilute alloy is 



V/Vi 



D 



1 + V/Vd ' 



(34) 



where diffusive speed Vd is a characteristic trapping velocity. For V/Vd 3> 1 
the segregation coefficient can be approximated by 



= l _ (l _ ke ) 



Vd 
V ' 



A comparison of the equations 1(33(1 and ((35(1 gives 



V D = 



D 



Al A {l - fc e ) RT 



df B df A 



(35) 



(36) 



b=l/2,T 



According to ((36(1 . the diffusive speed Vd depends on the diffusion coefficient, 
on the temperature, and on the difference of free energy densities of the pure 
materials. 



6 Numerical Calculations 

To validate the obtained solution, we will compare it with numerical solutions 
of the steady-state governing equations (JZJ) and JHJ. The material parameters 
used in the numerical calculations are given in Table 1. These parameters are 
similar to parameters employed in |S] and correspond to the alloy Ni-Cu [2] ■ 
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Table 1. Thermophysical properties used in calculations 



Parameter 


Value 




Reference 


La 


2350 


J/cm 2 


u 


Lb 


1725 


J/cm 2 


6 


T A 


1728 


K 


A 


T b 


1358 


K 


A 


D 


10" 5 


cm 2 / 


6 


a 2. 


8 x icr 5 


J/cm 2 


A 


Ia 6.< 


t8 x 1(T 8 


cm 


A 


V m 


7.4 


cm 3 /mol 





/" 


24 


cm/ (s-K) 


DDI 



The value of kinetic coefficient was set to \i — 24 cm/(s-K) according 
to the work JHj in which good agreements between the model for dendritic 
growth and experimental data have been achieved. The far-field concentration 
was set to Coo = 0.1. The parameter e is related to the surface energy a and 
to the interface thickness I a by equation e 2 = 6ctZa • The parameters Wa 
and Wb were chosen to be equal and are given as Wa — Wb — W — 12a /I a 
® 

Employing (J3J and J3J , a finite-difference approximation of the equations 
10 and (JHl on a uniform grid gives the following set of equations 



At Mie 2Z\ir Z\x 2 

TT (dg\ k _ W(c(3 B + (1 - c)/^ ) ( dp" 
6 



(37) 



„fe+i 



At 



V i 
M^7 



2Ac 



-i+l 



-24- 

Ar 2 



2Ax 



(38) 



g? = cf(i-c?) 

which determine the discrete functions of the solute and of the phase field 
at a point z = iAx on a (H l)th step of iteration with a time step At. To 
equalize the velocity V of the moving frame with the velocity of the interface, 
we introduce the equation 



dfB 

dcj> 



2Ax 



(39) 



(40) 
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where V s (r) is the shift velocity of the interface relative to the moving frame, 
Tf is an accessory parameter, r is a time variable corresponding to the it- 
eration process, and the interface coordinate zt is given by the condition 
cf>(zf) — 1/2. A finite-difference approximation of equation (|4U[) gives 



V 



fe+i 



V k 



eMi Zf 



At 



At 



(41) 



Results of numerical calculations of the equations H37 |) -l|41 |) are shown in 
the following figures. Fig. 1 shows a good agreement between the analytical 
solution [equations l|ll|) . I|2U|). (|23l) . and 124|) ] and the numerical results in 
the wide range of undercoolings. 



500 




15 20 

Undercooling AT, K 



Fig. 1. Kinetic curve "interface velocity V — initial undercooling AT" . The 
solid curve shows the approximated solution, data points denote the result 
of the numerical calculations. 



Fig. 2 shows the solute and the phase fields. At the interface, the solute 
field is smooth as against sharp-interface models, which gives a jump in con- 
centration. At low undercoolings, the analytical solution gives some excessive 
value of solute concentration in the vicinities of the interface. The series of 
numerical calculations given in Fig. 2 demonstrate the degeneration of the 
solute to uniform field with increasing interface velocity and undercooling. 
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Such behavior of the solute field is in accordance with the equation l|32Jl , in 
which the deviation from the far-field concentration is proportional to 1/V. 




Fig. 2. Solute profiles and the phase field for different values of the initial 
undercooling. The solid curve shows the approximated solution, data points 
denote the result of numerical calculations. 



7 Conclusions 

Planar solidification from undercooled melt was investigated by numerical 
and analytical methods in the limit of small impurity concentration. Com- 
paring the obtained results with those of Jf2 we note that while equilibrium 
partition coefficient (12(it has the same form as in work the expression 
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1)30(1 for the liquidus slope differs from the one given in the work 6 . However 
the evaluated equilibrium liquidus slope has very close values, m e — 306.2 K 
from equation ((30(1 and m e = 306.9 K from equation (A13) in 0. The dif- 
fusive speed has a weak dependence on the temperature, Eq. 1(36(1 . and can 
be accepted to be constant as well as the diffusive speed obtained in [H] . The 
numerical calculations presented here show a good agreement with the ana- 
lytical solution that indicate a sufficient accuracy of the expansion 1(10(1 - 1(12(1 . 

The obtained solutions simultaneously are valid for the bulk phases and 
the interface. Therefore they can be suitable for the purpose of numerical sim- 
ulations: for instance to construct initial splitting of adaptive finite-element 
grid and as a first approximation for the numerical solution of transcendental 
equations. 
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